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SUMMARY 


Volume-energy relations that can be determined from an equation of state are used 
to group many sets of heat-transfer data, ranging from data for liquids to data for gases 
and including the near-critical region. The volume - Gibbs energy parameter (or 
Helmholtz energy parameter) grouped these data better than did such other parameters 
as enthalpy, temperature, or internal energy. 


INTRODUCTION 


In relating volume-energy, or density-energy, parameters to turbulent flows 
(ref. 1), consider the conservation equations with fluctuations (refs. 2 and 3). The tur- 
bulent momentum and energy fluxes for steady flow in a tube become 

-T^ = p uV + V p'u' + p'uV ( 1 ) 

-q^ = p v'h' + V p'h' + pVh' (2) 

where u',v' are the velocity fluctuations and p',h’ are the density and enthalpy fluc- 
tuations, respectively. Symbols are defined in appendix A. The density, or specific 
volume, can be related to the enthalpy through an equation of state; 

p = p(h, p, xp (3) 

and the local density can be represented in terms of the average density by a Taylor 
series 


= p + -^Ah + -^AP + -^Ax. + . 
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or by using the volumetric parameters and fluctuating components (ref. 1) 

^^-/3„h' + kP' - f.x! 

^ ^ 11 


(4) 
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For single-component fluids at nearly constant pressure, the density fluctuations can be 
written as 


p' ^ 


(5b) 



The authors of 


where h' is the fluctuating turbulent energy ano (-9 In P/9h)p^^. 

reference 1 used these terms and the parameters Nup and or ip 2 group sev- 
eral heat-transfer data sets, ranging from low-temperature-liquid to high-temperature- 
gas data and including the near-critical region and high-pressure fluid states, where 
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i^'2 = 1 + ^2 ^w “ 
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and is of the McAdams type. Herein 


NUoaio = «-'523Re»-8pr“-'‘ 
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where properties are evaluated at the bulk temperature. The use of the Petukhov rela- 
tion is equally valid; 
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The equation of state representing the PVT surface can also be written several dif- 
ferent ways, such as 


P = p(<Pk» P. xp 


( 12 ) 


where 



^ = A = E - TS 

O 

= G = h - TS 


(Pg = E=h-PV ■] 



( 13 ) 



Does it follow that ip^, where (3^ =--- -(d In p/d(f^ , will also group these data? In ref- 
erence 1 it was shown that for k^l,2 the answer- is yes; and in this report we investi- 
gate k > 2, with the thought that, for some processes, it may be more convenient to use 
energy potential functions rather than, for example, temperature. 


THEORETICAL RELATIONS 

From an equation of state, such as equation (12), the density can be written to a 
first approximation as 


Ap ~ ^ Acp, -f ^ AP + Ax. (14) 

d(p^ 9P 9x. ‘ 

(no sum on k) where the Helmholtz, Gibbs, and internal energies were defined in the 
preceding equations/ 

In turbulent flow, since p p + p\ equation (14) represents the turbulent density 
fluctuation. For a process with a single-component fluid at nearly constant pressure. 


(no sum on k) where - -(9 In p/dcp^ and represents the turbulent energy 

fluctuation. The particular form to be used would depend on the variables of the energy 
equation. 

Thus we seek to determine the generalized parameter 


ipk - I + - ^k,h> 

(no sum on k) that will group heat-transfer data sets for k > 2 as was done for k < 2. 
We will need the following relations, which can be developed by using the Bridgeman 
tables (ref. 4); for example, at constant pressure 


9A _ / s ^ p^y. 9G _ = _ PV; 

9 1nV y 9 1nV 9 In V ’ 9 In V ( 3 ^ 


( 17 ) 
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for the simple fluids considered herein. 



RESULTS AND DISCUSSION 


The figures presented herein are based on the data sets delineated in reference 1 
and used to establish that Nu^ and or tp 2 can be used to group heat-transfer data. 

Gas -State Data 

For gases, to the first order is inversely proportional to temperature ^L.e., 

~ T~^^ . In reference 1, it was shown that for k < 2 

Nu^ = C4>~^ (18) 

where 0.95 < C < 1.1 and 0.45 < n < 0.6. Representative values were given for C 
as 1 and for n as 0.55 for fluids and 0.50 for gases (ref. 1). As 4> approaches 1, 
the correction to the conventional McAdams equation becomes insignificant. (Note that 
for a gas, T^/T^^.) 

In figure 1, Nu^^ is shown as a function of where k - 1, 5, for the data sets of 
references 5 and 10, which represent four gases - helium, hydrogen, air, and carbon 
dioxide. The grouping appears to follow equation (18). The figure presents 1345 data 
points that appear to be represented by C~1.05 and n~0.5 for k = l,2,5 and n~0.4 
for k = 3,4. The figure is intended to show that these data are grouped by (k = 1,5) 
and not to develop specific values of C and n. The individual data sets are discussed 
in appendix B for k = 1. 


Fluid-State Data 

In the fluid state, which includes the liquid and dense-gas regimes, the thermo- 
physical properties vary significantly and the situation is more complex. The data for 
fluid hydrogen are plotted in figure 2 for low-pressure hydrogen flows through a 0.478- 
centimeter- (0.188-in.-) diameter tube, where buoyancy effects were considered negli- 
gible (ref. 10). The high-pressure data sets of reference 11 are plotted in figure 3. 
The high-pressure data tend to follow equation (18); for the low-pressure data, tp^ and 
ip^ more closely follow equation (18), but ip^ appears to be cubic in log - log 4> 
for ip > 1.2. 

Fluid-nitrogen data (ref. 12) are plotted in figure 4. Fluid-oxygen data (refs. 13 
and 14) for low pressure and high pressure are plotted in figures 5 and 6, respectively. 
The nitrogen and low-pressure oxygen data tend to follow equation (14) , with the nitro- 
gen data 30 percent low. The high-pressure oxygen data appear to be twice the values 
predicted by using equation (18) at low ip. For both these data sets, it is felt that the 



inside surface temperatures have not been properly corrected. Correction could result 
in a 25- to 30-percent increase in the nitrogen data and a lowering of the oxygen data. 

Figure 7 (ref. 15) represents helium data for k = 1 . No other comparisons were 
made. 

Although each groups these data sets, the Gibbs and Helmholtz jp^ func- 
tions group the data better than temperature jp^, enthalpy ip^, and internal energy ip^ 
in that order. This judgment is based on a comparative analysis for many plots; over- 
lays of the various data sets; an analysis of data scatter; the fluid state or regime 
where this scatter occurred; and the relative steepness of the slope of the data group- 
ing, all of which can contribute heavily to predictive errors. Such a qualitative trend 
suggests that the turbulent energy equations should perhaps be developed in terms of the 
Gibbs or Helmholtz functions, depending on the process involved, for a better under- 
standing of the dynamics. 

The hydrogen, oxygen, and helium data of figures 2, 6, and 7, respectively, indi- 
cate an enhanced heat transfer at low more notably in the critical region. An in- 
crease in the deviations at high ip^ is noted for hydrogen ^ig. 2) . The divergence of 
Nu^ at high ipy^ may be caused by flow disturbances (oscillations) for at high ip^^ the 
flow is quite sensitive to small disturbances and oscillations are commonplace (e.g., 
ref. 10). Therefore, it may become necessary to incorporate /c = 9 In p/9P into 
that is, 

i/'k = 1 ^ Acpk + /c AP (19) 

These results again indicate that data in the near-critical region can be grouped 
better by using jp^ or !p^ and that high-pressure and gas-phase data can be satisfac- 
torily grouped by using 


2 

For the data herein, L/D was not considered important; however, where L/D is 
important, two forms suggested in ref. 1 can be applied: 

Nu^ = 
or 

Nuj^[[l+ (2D/L)]j =ip^^ 

where 


m = n - (1 + n) D/L. 
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SUMMAHY OF RESULTS 


An analysis of the turbulent transport fluxes, and heat-transfer data from several 
sources reveal the following results: 

1. Volume-energy, or density-energy, parameters such as ^-9 In p/dcp^ that 

are readily obtained from an equation of state can be related to turbulent flow density 
fluctuations, where p is density and is an equation-of-state variable. 

2. Volume-energy parameters can be used to group heat- transfer data sets ranging 
from liquid-state to gas-state data, including the near-critical region. The Gibbs or 
Helmholtz functions appear to group data better than other thermodynamic parameters; 
however, outside the near-critical region, volumetric expansion (k = 1) can be used. 

3. These results show that some insights might be gained by reformulating and 
solving the turbulent flow equations in terms of the Helmholtz or Gibbs functions . 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 31, 1979, 

506-16. 
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APPENDIX A 
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SYMBOLS 


Helmholtz energy 
constant 

specific heat at constant pressure 

diameter 

internal energy 

Gibbs energy 

enthalpy 

length 

exponents 

Nusselt number 


pressure 
Prandtl number 
turbulent energy flux 
Reynolds number 
entropy 
temperature 

axial velocity component 

specific volume 

radial velocity component 


fraction of species i 

isoparametric volumetric expansion, - (9 In p/dcp.) 

>k 

isoparametric composition coefficient, -(9 In p/9x j 
isoparametric compressibility, (9 In p/9P)' 




defined by eq. (11) 
chemical potential 
defined by eq. (11) 
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p density 

turbulent momentum flux 

cp equation-of-state or fundamental-equation variable 

ijj defined by eq. (16) 

Subscripts: 

b evaluated at bulk conditions 

k index 

R reduced, see eq. (6) 

w evaluated at wall conditions 

Superscripts: 

( ) average 

( fluctuating 
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APPENDIX B 


GAS-STATE DATA SETS 

The data of individual investigators that constitute figure 1 require further analysis. 
These data sets appear to be self-consistent but exhibit significant deviations when com- 
pared by using such common parameters as those of equation (18) (i.e, , Nu^^^ and ip^). 


Gaseous Hydrogen 

McCarthy (ref. 5) acquired an extensive data set for gaseous hydrogen flowing in 
tubes over a wide range of the control parameters, namely, the inlet conditions, mass 
flow, heat flux, and wall temperatures. The high thermal driving potential (T^ - T^) 
was in some cases achieved by precooling the gas before running it through the Joule- 
heated test section. Wall temperatures to 1260 K were achieved. These data, plotted 
as figure 8, are characterized by equation (18) where n = 0.5 and C = 1.1. 

Weiland (ref. 7) also took an extensive set of data, although the range of heat fluxes 
was not as extensive as McCarthy's (ref. 5). The test-section L/D ranged to 237. 
Weiland also used a precooler to condition the gas. The axial temperature profiles 
showed significant departures from the conventional monotone turbulent-flow profile. 

The data are plotted as figure 9. For these data, n == 0.6 and C - 1.0. 

Taylor (ref. 6) took an extensive set of data with very high wall temperatures (to 
3000 K) by using a Joule-heated molybdenum tube. Under these conditions many ques- 
tions as to accuracy in measurements and thermophysical properties enter. And even 
though no abnormalities or dissociation were reported, still some question exists at the 
higher values of These data are plotted as figure 10. Here the data follow one 

form of equation (18), namely, n = 0.5 and C = 1.0 for < 5, and appear to have 
little influence on ip^ for values greater than 5. 

Hendricks, et al. (ref. 10) acquired a limited set of data for five tube diameters, 
ranging from 0.478 to 1.28 centimeters (0.’l88 to 0,504 in.). In general the driving po- 
tential is low, and the data should not be given the same weight as the data of references 
5 to 7. These data are plotted as figure 11. 


Gaseous Helium 


To investigate the effects of an inert gas flowing in a Joule-heated tube, McCarthy 
(ref. 5), Weiland (ref. 7), and Taylor (ref. 6) also took extensive sets of helium data in 
their respective apparatus. These data are plotted as figures 12, 13, and 14, respec- 


9 



tively. The consistency of the helium and hydrogen data from each investigator is ap- 
parent and would appear to rule out dissociation effects in the data of Taylor (ref. 6). 

It should also be pointed out that McCarthy (ref, 5) and Weiland (ref. 7) used log- 
average values and Taylor (ref, 6) used local values, i^lso suspect are the thermophys- 
ical properties routines, as many temperatures are well beyond the approved range 
(ref. 4). 


Other Gases 

Since the properties of gaseous air are quite similar to those of nitrogen, the prop- 
erties of nitrogen were used to reduce the set of air data taken by Wolf (ref. 8) . These 
data, although grouped by equation (18), appear to have some geometry effects not noted 
in the other data sets. These data are plotted as figure 15. Wolf (ref. 8) also took data 
for gaseous carbon dioxide. These data are plotted as figure 16, These results also 
show trends with geometry and may indicate these data to be apparatus dependent. 

The gaseous air data set of Humble (ref. 9) was also reduced by using the proper- 
ties of nitrogen. These data are plotted as figure 17 and show more scatter than the 
previous sets. 

In general the gaseous data are grouped quite well by equation (18), an indication 
of the usefulness of the density-fluctuation terms (eq. (15)). And for gases 0.95 < C 
<1.1 and 0.5^n<0.6, 
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(c) Helmholtz energy. 


Nusselt number for fluid oxygen (ref. 11) as a function of volume-energy parameters. 
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Figure?. - Reduced Nusselt number as a function of iji ^ for fluid 
helium (ref. 13). 





Figure 8. - Reduced Nusselt number as a function of the 
volumetric expansion parameter il/i for gaseous state - 
hydrogen data from reference 5. 
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Figure 9. - Reduced Nusselt number 
as a function of the volumetric 
expansion parameter for gas- 
eous state - hydrogen data from 
reference 7. 
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Figure 10. - Reduced Nusselt number as a function 
of the volumetric expansion parameter for 
gaseous state - hydrogen data from reference 6. 
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Figure 11. - Reduced Nusselt 
number as a function of the 
volumetric expansion param- 
eter iJ/Y for gaseous state - 
hydrogen data from refer- 
ence 10. 
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Figure 12. - Reduced Nusselt number as a function of the 
volumetric expansion parameter 0i for gaseous state - 
helium data from reference 5. 
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Figure 13. - Reduced Nusselt number 
as a function of the volumetric 
expansion parameter i//-^ for gas- 
eous state - helium data from ref- 
erence 7. 
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Figure 15. - Reduced Nusselt number 
as a function of the volumetric ex- 
pansion parameter ipi for gaseous 
state - air data from reference 8. 


Figure 14. - Reduced Nusselt number as a function of the 
volumetric expansion parameter for gaseous state - 
helium data from reference 6. 
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Figure 16. - Reduced Nusselt number 
as a function of the volumetric ex- 
pansion parameter 11/^ for gaseous 
state - carbon dioxide data from 
reference 8. 
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Figure 17. - Reduced Nusselt number 
as a function of the volumetric ex- 
pansion parameter for gaseous 
state - air data from reference 9. 
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